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Phase Shift Effects in Fabry-Perot Interferometry 
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(January II, L960) 

A method is demonstrated for utilizing in Fabry-Perol interferometry the data on 
reflection phase shift dispersion obtained from fringes of equal chromatic order. Unknown 
wavelengths can be calculated from the Fabry-Perol patterns obtained with a large etalon 
spacing, even without prior knowledge of the phase shift of the reflecting surfaces. When 
the theoretical phase shift as a function of wavelength is known approximately, then the 
correct orders of interference can be determined for both the Fabry-Perot fringes and fringes 
of equal chromatic order. From the wavelengths of the latter the phase shift dispersion can 
be measured to an accuracy of about 10 A. The method is especially useful for reflectors 
with large dispersion of phase shift, such as multilayers. Results in the visible spectrum 
are reported for aluminum films and a pair of dielectric 15-layer broadband reflectors. 



1. Introduction 

In the field of Lnterferometric Length and wave- 
length measurement it has long been recognized that 
the dispersion of phase shift on reflection must be 
considered. In Fabry-Perot interferometry, pro- 
cedures for determining and/or correcting for phase 
shift dispersion have been described by Fabry and 
Buisson [I], 2 Eversheim [2], Meggers [3], Bauer [4], 
Jackson [5], Meissner [6], Barrell and Teasdale- 
Buckell [7], and by Rank and Bennett [8]. 

In wavelength intercomparison, the* phase shift 
dispersion enters the calculations at two points: 

(1) The determination of the etalon spacing, and 

(2) the calculation of the unknown wavelength. 
The purpose of this paper is to describe how the 

accurate measurement of phase shift dispersion by 
fringes of equal chromatic order (feco) [9, 10, 11] can 
be used with Fabry-Perot interferometry in the 
determination of the etalon spacing and in precise 
measurement of unknown wavelengths. 

1.1. Determination of Etalon Spacing 

The method of exact fractions [12, 13] is generally 
used to determine the etalon spacing from the 
observed interference pattern. For three known 
wavelengths, X a , X 6 , X c , the fractional orders J a , f b , 
and f c , are calculated [6, 14] from the diameters of 
the circular fringes. Then the integral order num- 
bers, m a , m b , m c , and the exact etalon spacing, /, 
are found by trying many different integers for 
m a , m h) and m c until the products (m a +/ a )\ a , 
(m 6 + / 6 )X 6 , and (m c +f c )X c are all essentially equal. 
The etalon spacing is then set equal to 

K(m a +/ a )X a =K(m 6 +y 6 )X 6 =K(^c+/c)Xc (1) 

This method obviously ignores the phase shift. 
The quantity }^(m+/)X really represents the optical 
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spacing of the etalon plates, which is the physical 
separation plus the phase shift on reflection for 
the particular wavelength. Since in general the 
phase shift is a function of wavelength, the quantity 
){{m-\-f)\ is not constant as implied in eq (1) but is 
a function of wavelength. It is fortunate that for 
the metals aluminum and silver which are commonly 
used on Fabry-Perot mirrors, the phase shift disper- 
sion is such that it is possible to find integers m aj ra&, 
and m c which satisfy eq (1). When two different 
etalon spacings are determined by this method 
correct wavelengths values can be calculated. 

Multilayer reflectors [15, Hi, 17, 18] are finding 
considerable application in interferometry because 
of their low absorbanee and high reflectance. For 
the common type of high reflection coating consisting 
of alternate high and low index layers, (Mich of 
quarter-wave thickness for wavelength X, the phase 
shift dispersion is small in the wavelength region 
near X. Outside this spectral region the phase 
shift dispersion is much larger than for silver and 
aluminum mirrors. A 15-layer broadband high 
reflectance multilayer described by Baumeister and 
Stone [19] has a very large phase shift dispersion 
[20]. In such cases of high phase shift dispersion 
it is possible that integers cannot be found which 
will satisfy eq (1), or wrong integers may be selected. 

Rank and Bennett [8] have described a method 
for using calculated values of the phase shift to 
correct the order number of interference in a Fabry- 
Perot etalon. Baumeister and Jenkins [20] have 
described another method for employing calculated 
phase shifts and have applied it to a nine-layer 
coating. The latter authors have also made measure- 
ments of the phase shift produced by the 15-layer 
broadband multilayer using Fizeau fringes [1]. The 
precision of this method was limited by the broad 
two-beam fringes used for comparison. 

1.2. Calculation of the Unknown Wavelength 

For work of highest accuracy, the above method 
of exact fractions is used only to determine the-- 
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correct order numbers. The e talon spacing and 
the unknown wavelength are computed from a 
standard wavelength X s as follows. The method 
eliminates phase shift dispersion errors. 

The fractional orders of interference, j x and f s are 
obtained for the unknown wavelength, X^, and for 
the standard wavelength, which can be one of the 
wavelengths used to determine the etalon spacing. 
The equations for interference at normal incidence 
are 



and 



(m x +f x )\ x =2t+2<K\ x ) 

(m s +f s )\ s =2t+2d>(\ s ). 



(2) 
(3) 



Here <£(X) is the phase shift in units of length; it is 
the additional optical path introduced on reflection. 
There are as many as four unknowns in these 
equations: X^, t, 4>(\ x ), and <£(X S ). (The unknown 
wavelength, X^, must be known in advance with 
sufficient accuracy to determine the order number, 
m x .) 

If another etalon spacing, t', is used, two more 
equations are obtained 



(m' x +f x )\ x =2t'+24>(\ x ), 
(7ns + r s )\=2t'+2$(\). 



(4) 
(5) 



The order number ml can be determined by the 
method of exact fractions described above. Again 
m x is determined from the approximate value of X^. 
The four equations can be solved [6] for \ x by forming 
the difference between the two equations involving 
X x , also between the two equations involving X s , and 
equating the differences. 



X-r — X, 



nis+fs 



-m. 



-/: 



m x +f x —m x —f x 



(6) 



The phase shift terms are thereby eliminated. The 
difference in etalon spacing, t— t', is obtained by 
subtracting eq (5) from eq (3) 



t — t' = %\ 8 (m 8 +f 8 - 



■m 8 



-/0. 



Variations of the above method of analysis have 
been used by Meggers [3], Jackson [5], and Meissner 
[6]. The variations make analysis of an extended 
series of data more convenient, but fundamentally 
they depend on the measured quantities, f s , /J, f x , 
and /i in the manner shown in eq (6). 

The order numbers m s , m 8y m x , and m x are integers 
which can be determined without ambiguity in cases 
where <t>(\ 8 ) and 4>(k x ) can be adequately approxi- 
mated. The accuracy with which X^ can be calcu- 
lated from eq (6) therefore depends on the accuracy 
with which the fractions, /, can be measured. In 
particular, the error in X^ is given approximately by 



J */. Y+ ( *r. v + ( a/, y, / a/; y 

\m s —m'J \m s —rn' s J \m x —m x / \m x —m' x ) 

(7) 

where AX^, Af s , Aj s , Af x , and Af x represent the stand- 
ard deviations in the values of X^, f s , f SJ f x , and j x , 
respectively. 

If the phase shift dispersion can be determined by 
independent measurements, then the unknown wave- 
length can be calculated from eqs (2) and (3) 









x*= 



(m,+/,)A f +20(A I )-2*(A,) 



m x +f x 
The error in A^ is given approximately by 



(8) 






m 



E 1 =('w&ve incident on surface 2) 

=A 1 e i{<aT - KX *\ 
E 2 = (reflected wave) 



_/A/A 2 /A/A 2 / 2A4>(K) \ 2 ■ / 2A»(X,) V (g) 
\m s J \m x ) \ m s \ s ) \ m s \ s ) 

Comparison of eqs (7) and (9) shows under what 
conditions the second method (eq 8) gives more 
accurate results than the first method (eq 6). For a 
given value of A/ s , the value of Af s /m s in eq (9) will 
be smaller than Af s /(m s -m s ) in eq (7). However, 
the difference need not be great if m s is small com- . 
pared to m, s (that is, if the smaller etalon spacing, 
t', is small compared to the large etalon spacing, t). 
The greatest gain is realized if 2A<t>(\ s )/m s \ s can be 
made significantly smaller than Aj 8 j{m s — m s ), simi- 
larly for the quantities involving the unknown wave- 
length, \ x . Since it is possible to have m s —m f 8 ^m 8 , 
this requirement can be written as 2Ac/>/X<A/ / . 
The limit on determination of the fraction/' from the 
diameters of the circular fringes is about Af = 0.01 
[21,22]. It will be shown that when fringes of equal 
chromatic order are used to determine the phase 
shift dispersion, the quantity 2A</>/X can be con- 
siderably less than 0.01. 

2. Theoretical Calculation of Phase Shift on 
Reflection 

The convention used in this paper is that the phase 
shift on reflection, <f>, represents an increase in optical 
path. The angular phase shift, e = 27r<£/X, enters 
the expressions for the reflected waves as follows 
(6=0; see fig. 1) 



(10) 



(ID 
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E 3 =(E 2 iucident on surface .1 from inside) 

= A 2 e i{ " T - Kr *- € - d) , (12) 

Z?4= (reflected component of E$) 

= A 4 e H " T - KX >-*>- e i- Kt \ (13) 

E & = (wave E 4 incident on surface 2) 

=A A e i{u}T - K ^-^-^- 2Kt \ (14) 

Here a> is the angular frequency, T the time, k = 2t/\. 




Figure 1. Schematic diagram of interferometer. 
The reflecting surfaces are denoted by 1 and 2. For feco and for the center of the 
Fabry-Perot Interference pattern the angle of incidence, 0, is zero. 



The condition for constructive interference in 
transmission is that the optical path difference be- 
tween two successive rays should be an integral 
number of wavelengths. 

2irn=avgE l — &Yg E 5 =coT— kx 2 — (coT— kx 2 — e 2 — ci — 2d) 

= e 2 J re\ J c2Kt 

n\=cl) 2 +(t) 1 + 2t. (15) 

From eqs (10) and (11) we have 

E l A, e 

For metals or for dielectric multilayers, theory can 
be used to calculate the quantity 

— = — - p l h. 

E l A 

Therefore the relation between the desired angular 



phase shift, e, and the calculated value, 8, is 



€ = J>27T — 8 



(16) 



where v is an integer. 

For a metal reflector the angular phase shift, 8, 
is calculated from the index of the incident medium, 
fjL , and the optical constants, p and /', of the metal. 
The expression for the complex amplitude reflectance 
is [2U] 



pe t{ 



_ Ho—ii+ik _ 



j4—n 2 —k 2 +i2iJL k 

Mo + M 2 + ^ 2 + 2moM 



(17) 



The tangent of the phase shift 8 is then found by 
taking the ratio of imaginary part to real part: 



tan 8= 



2 nJc 



2 -¥ 



(18) 



The quadrant for 8 is determined by taking into 
account the signs of the real and imaginary parts of 
eq (17). Thus for /z 2 +£ 2 >/io (the case for all metals 
commonly employed in interferometry, when the 
separation medium is air or vacuum) the real part is 
negative and the imaginary part is positive, which 
means that 8 is in the 4 second quadrant. The phase 
shift = eX/27r=(v27r — 8)\/2tt therefore represents 
an increase of optical path of between A/2 and 3X/4. 3 
Figure 2 shows the phase shift for aluminum calcu- 
lated from the optical constants [24, 25]. 

For a multilayer the phase shift can be calculated 
from electromagnetic theory by one of several meth- 
ods. An accurate method employs multiplication 
of matrices, each matrix representing a single layer 
[26, 27]. An exhaustive treatment of the subject 
has been given by Abeles [28], but the outline given 
by Koehler [29] is adequate for many problems. 

206-I0 3 
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"i — 1 — 1 — 1 — 1 — 1 — r 
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Figure 2. Dispersion of phase shift on reflection for aluminum 
films. 
The theoretical values were calculated From the optical constants as given by 
Schulz [24, 25]. For the upper experimental curve the data from table 1 and the 
Integer r/=31 were used. Since this curve is parallel to the theoretical curve, it 
follows that f/ = 31 is the correct order number for the fringe at X =6593.7A. For 
the lower experimental curve the integer q='3Q was used . 



3 The integer v=l was selected for simplicity. Some other authors select i>=0, 
which makes represent a decrease of optical path of between X/4 and X/2. For 
an infinite wave train the two interpretations are physically identical. They do 
lead to different order numbers of interference. 
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In dealing with phase shifts on reflection from 
multilayers the concept of nodes is very useful [20]. 
For light reflected from a plane interface between air 
and a high index dielectric, the phase shift is <f> = \/2, 
and there is a node at the surface. For </>>\/2, as 
with a metal, the node is inside the metal. For 
<f><^\/2, the node is on the air side of the surface. 
With multilayers of the common type, having a high 
index layer on the outside, a node for the tuned wave- 
length, X, falls at the outer surface. For the broad- 
band multilayer it has been shown [20] that there is 
one wavelength for which a node falls at the outer 
surface. At longer wavelengths the node is outside 
of this surface, and vice versa. Figure 3 shows the 
dependence of phase shift on wavelength. 




- 334-I0 3 
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Figure 3. Dispersion of phase shift on reflection for 15 layer 
broadband reflecting films. 

The angular phase shift, S, as given hy Baumeister and Jenkins [20] was con- 
verted to linear phase shift </> by means of eq (16) in which v=\ and the relation 
<t> = cX/2tt. The experimental curve has the same general shape except that it is 
shifted to longer wavelengths. Wavelengths from table 2 were used. 



3. Experimental Procedure 

Data on the phase shift dispersion were obtained 
by means of fringes of equal chromatic order [9, 10, 
11, 30]. As shown in figure 4 collimated white light 
illuminated the interferometer. The very small 
spacing, 10 to 30 u, was accomplished by means of 
aluminum foil. The transmitted light entered a 
spectrograph, with the image of the interferometer 
surface being focused on the slit. As illustrated in 
figure 5, bright fringes appeared in the focal plane 
of the spectrograph for those wavelengths at which 
constructive interference occurred. After photo- 
graphing the fringes along with an iron arc reference 
spectrum, the wavelengths of the fringes were meas- 
ured. The method of calculating the phase shift 
dispersion is given below. 

Fabry-Perot interference patterns were obtained 
using invar spacers of 2 and 5 mm, also fused quartz 
spacers of 20-mm length. The light source was a 
natural krypton lamp, krypton being selected be- 
cause of its relatively sharp and well-distributed 



lines. The etalons were used in air maintained near 
20° C and 10-mm water vapor pressure, and correc- 
tions to the wavelengths were made by the method 
described by Bruce [31]. 

From the photograph of the interference pattern, 
the ring diameters were measured for several wave- 
lengths. For the 2-mm spacer, two rings were meas- 
ured at each wavelength; for the others, four rings 
were measured. The fractional orders, /, of inter- 
ference at the center of the patterns were calculated 
with a least squares procedure described by Meissner 
[6]. 
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nXo=2^+20(Xo), 

(n+l)\ 1 = 2t'+2<l>(\ 1 ), etc., 



(19) 
(20) 



X and X, are the wavelengths of adjacent fringes, 
\ >Ai. The order number for the fringe at X is n, 
and for each successive fringe the order number is in- 
creased by one. t' is the physical separation of the 
reflecting surfaces, and <£(X) is the phase shift in 
units of length at each surface. 

In these equations the left sides can be evaluated 
only at the wavelengths X , X b etc., at which there are 
fringes. However the quantity on the right side is a 
continuous function of X which is conveniently de- 
noted as r n (\), 



Tn (\)=2t'+2<t>(\). 



(21) 



The function r n (\) is then evaluated by plotting 
n\ versus X , (n-\-l) At versus \ x , etc., and connecting 
the points by a smooth curve. This curve has the 
same shape as the 2</>(X) curve but is displaced verti- 
cally by 2t'. 

With the Fabry-Perot etalon the interference equa- 



Figure 4. Interferometer for fringes of equal chromatic order. 

A. White light source. B. pinhole, C. collimator, Y>. beam-divider (for reflected 
fringes), E. flats with metal or multilayer reflecting films, F. achromatic doublet, 
(J. interferometer base, and II. spectrograph. 

For Fabry-Perot Interferometry the setup is changed as follows: The line 
source is placed at the focal point of the collimator, C. A larger spacer is intro- 
duced between the flats, F. A lens with longer focal length than F is used to 
focus the fringe pattern on the spectrograph slit, which in this case is wide. 



4. Equations for Fringes of Equal Chromatic 
Order and Fabry-Perot Fringes 

For the wavelengths X , X ]? etc., at which feco occur 
in the spectrograph, the normal incidence interference 
eq (15) can be written as 
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FIGURE 5. Fringes of equal chromatic order with iron reference spectra (positive print). 
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a. Aluminum films, reflection fringes, 10-/1 spacer; 

b. aluminum Alms, transmission fringes, LO-/* spacer; 

c. 15 layer broadband films, transmission fringes, L 7-/i spacer. 

In C, the only obvious Indication Of the large dispersion of phase shift is the f id tint t he fringe spacing docs not decrease monotonically with increasing wavelength. 
In the red the last few fringes are more widely spaced than those Immediately preceding. Where several fringes are overexposed in the red the standard deviation 
in determination of wavelengths is about 0.4 A. whereas in other portions of the spectrum it is from 0.1 to 0.3 A. 



tions for the centers of the ring pattern a I known 
wavelengths X a , X 6 , X r , are: 



(ma+fa)K=2t + 2<t>(\ a ), 

(m b +f b )\ h =2t+2<t>(\ b ), 
(m c +f c )\ c =2t+2<t>(\ c ). 



(22) 
(23) 

(24) 



Here again the quantity on the right side is a continu- 
ous function of X, the eurve being displaced upward 
from the curve r n {\) by the amount 2(t— t'). There- 
fore, in general the points (ra fl +/ a )\ tt , (m 6 +/ 6 )X 6 , 
and (m c + / c )X c lie along a curve w T hich is parallel to 
the r n curve. 

5. Methods of Analysis 

5.1. Method I 

The first method for obtaining the order numbers 
n and m is used in cases where the phase shift 
dispersion is known approximately, as for example 
with metallic films with known optical constants for 
various wavelengths. For the jeco an integer p is 
selected and the curve 1 t>(X) is determined by 
plotting /A versus X , (p J r \)\ versus Xi, etc. If it 



does not have the same shape as the known 2</>(X) 
curve, then other integers are tried until an integer, 
n, is found which makes r n (X) essentially parallel to 
2(/>(X). This procedure is illustrated in figure 2. 

Similarly for the Fabrv-Perot data, integers m ai 
m b , m c are found such that when {m a +ja)\ a is plotted 
against X a , (m 6 +/ 6 )X/, against \„, etc., the points fall 
along a curve parallel to r n (X). In order to measure 
an unknown wavelength, X T , in terms of a standard 
wavelength, X. s , a modification of eq (8) is employed. 
Although the absolute values of <p{\ x ) and <t>(\ s ) have 
not been determined, it can be seen from eq (21) 
that the difference, 2<f>(\ x )—2<t)(\ s ) i is given by 



2</>(X,)-2<KX,) = t w (X j )-t„(X 5 ). 



Therefore eq (8) becomes 



_ (ra g +/ a )X, + T w (Xg) — T n (\) 

m x +f x 



(25) 



(26) 



The method has been described in terms of 
graphical procedures. A possible variation utilizes 
numerical calculations similar to those in the usual 
method of exact fractions. For wavelength X„ the 
value of T n (\ a ) is found from the curve or by means 
of a numerical interpolation of the r n (X) data. An 
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integer r a is selected and the product (r a -f/ a )\ a is 
formed, also the difference (r a +/ a )X a — r n {\ a ). This is 
done for all the integers, r a , which lie within the range 
of possibility. Similarly different integers are tried 
for r b and r c until a set is found which satisfy the 
equation 

{r a +Ja)K— T n (K a ) = (r b +j b )\ b — T n (X 6 ) 

= Oc+jQX c — r w (X c ) 

Once the correct integers have been found by this 
technique eq (26) can be used to compute unknown 

x\m \rnl on erf h« 

5.2. Method II 



wavelengths. 



The second, more general, method can be used 
when the phase shift dispersion is unknown. An 
arbitrary integer, q, is selected and the curve r g (\) is 
determined as before. Then for the Fabry-Perot 
data, different integers are tried for r a , r b , and r c , 
until a set is found for which the points (r a +/ a )X„, 
(r b +f b )\ by and (r c + / c )X c fall along a curve parallel to 
T g (\). That is, 

(;r a +j'a)K—T Q (\a)=(r b +j b )\ b —T Q (\ b ) 

= (r c +f c )\-T a (\ e ) (27) 

In appendix 9.1 it is shown that when the integers 
rj are chosen in this manner then we have 



(rj+fj)h-T Q (h) = 2(t-t') 



(28) 



regardless of the value of q. Here the subscript j 
represents the subscript a, b, c, s, or x. 

The final step is to determine an unknown wave- 
length \ x in terms of a standard wavelength, X s , and 
the measured fractions j x and f s . The integers r s and 
r x are determined in the same manner as r G , r by 
and r c . When these quantities are substituted in eq 
(28) two equations result which can be solved for \ x : 



X*= 



z (rs+Js)K+r Q (\ x )-r q (\ s ) t 

r x +fz 



(29) 



Just as in the standard method, the wavelength 
\ x must be known in advance with sufficient accuracy 
that the proper integer, r XJ can be selected. The 
correct integer, r x , is the one which makes the 
quantity (r x +/ z )X x fall closest to the curve (r^-f^Xj 
versus \j. 

6. Results 

6.1. Evaporated Aluminum Films 

Table 1 gives the wavelengths of jeco obtained in 
transmission with a 10- fx spacer, the aluminum films 
having a transmittance of 5 to 8 percent. In prin- 
ciple the phase shift should be calculated taking into 
account the finite thickness of the aluminum films 
[32]. However, figure 2 shows that for the thicknesses 
used, the phase shift values calculated from the bulk 
constants yield a curve quite parallel to the experi- 



mental curve. Method I was used to determine the 
order number, n, as illustrated in figure 2. 

For the 2-, 5-, and 20-mm spacers the integers, m, 
were found as described in 5.1 and illustrated in figure 
6. With each spacer the points have the same verti- 
cal displacement from the curve r n , within experi- 
mental error. 

Table 1. Fabry-Perot and feco data obtained with aluminum 
films 

Ambient conditions for Fabry-Perot data, 20-mm, 5-mm spacers: 21.2°C, 754.6- 
mm Ilg barometric pressure, 10-mm vapor pressure; 2-mm spacer: 22.0°C, 758.8- 
mm, 10-mm vapor pressure. 
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As an example, the wavelengths 5570.34 and 
4502.39 are considered being known approximately 
to the nearest 0.02 A, this precision being necessary 
to establish the order numbers. The more precise 
values are to be determined by interferometric com- 
parison with the standard line at 6056 A. From 
figure 6, 

73i (5570) -T3i (6056) =203,350-203,830= -480 A. 

Therefore from eq (26) we have 

_ 66267.895X6056.1779-480 
x ~ 72047.649 

= 5570.3369 A. 

Similarly for X :r =4502 A 

66267.895X6056.1779-153 
x ~ 89136.737 

= 4502.3931 A. 

These values agree to experimental accuracy with the 
wavelengths for ambient conditions, as given in 
table 1. 

In illustrating method II, rather than choose an 
arbitrary integer, we choose an integer, q, such that 
r q is as nearly constant as possible. As shown in 
figure 2, t 30 has only a very slight negative slope. 
The integers, r h which make the points (rj+fj)^ lie 
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Figure C>. Optical path d /'])'< •retires: alum in inn films. 

A constant was subtracted from the (m,+/,)X,- values to bring them on the same 
scale as the tji curve obtained from transmission /«co. For the 20-, 5-, and 2-mm 
spacers the constants were 4011 1!MI) :| . 99698-10 8 , and 3969M0 8 A respectively. For 
the 5-mni spacer (open circles) several differenl sets of Integers, m,-, were tried 
together with the appropriate fractions, /,. One set Lives values of (TO,-+//)X* 
which lie along the dashe.f! line, which is parallel t()T3i. The other sets give points 
which obviously do not lie along lines parallel to T31. For the 2- and 20-mm 
spacers the correct integers plus fractions are given in table 1, and the optical path 
differencesjare plotted above. 

along a curve parallel to r 30 are one less than the 
corresponding integers m h in table 1. That is, for 
\=6456, r,=62,159. 

Then for the unknown line, X^ — 4502.39 A 

x 66266.89 X 56056.1779-197770+197,800 



89135.737 



=4502.3932 A. 



Thus with aluminum films if a certain integer, q, is 
chosen the r ff (X) curve is nearly constant and the 
phase shift correction is very small. Furthermore 
the quantities (^+/^)X^ are nearly constant. In the 
standard method of exact fractions the integers are 
selected so that the quantities (rj-\-fj)\j are equal. 
Therefore the method of exact fractions as normally 
used gives integers, r h which are each one less than 
the "true" order numbers, %. 4 This introduces no 
error into the determination of unknown wavelengths 
because only differences in order numbers for two 
different e talon spacings are used in the final 
calculation. 

6.2. Fifteen-Layer Broadband Reflectors 

Table 2 gives the wavelengths of feco obtained in 
transmission for a broadband multilayer [19] made 
from alternate layers of cerium dioxide and mag- 
nesium fluoride [33]. The measured reflectance was 
between 83 and 94 percent over the region 4500 to 
7000 A. In figure 3 the theoretical phase shift curve 
was obtained from values published by Baumeister 



* The word "true" is In quotation marks because some authors use another 
convention for sign of the phase shift, which results in another order number. 
See section 2. 



Table 2. Fabry-Perot and feco data obtained with 15-layer 
broadband ?nultilayers 



Amiuent conditions for Fabry-Perol data: 19.5° C, 757.05 mm Ilg barometric 


pressure, 6.5-mm vapor pressure. 
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4830. 4 
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5570. 3222 
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4693. 1 


71 
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4502.3809 


22185.150 


99,885,996 


4563. 4 


73 
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7107. (i() 


39, 925, 841 








4502. 3809 


SSI 17. 57 


39, 925, 178 








4273. 9952 


9341.38 


39, 925, 013 



and Jenkins [20]. The experimental phase shift 
curve confirms the general shape of the theoretical 
curve. The displacement of about 400 A to longer 
wavelengths indicates that the individual layers are 
somewhat thicker than t he* theoretical values [19]. 
This is confirmed by the location of the reflectance 
maximums and minimwms. 

For the 2-, 5-, and 20-mm spacers the integral 
orders of interference, m, were obtained as described 
in section 5.1. For each spacer, the total optical 
path differences (m<+/^)X^ plotted against X^ fell 
along a curve parallel to the r 48 curve, as shown in 
figure 7. 
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Figure 7. Optical path differences: 15 layer broadband films. 

T48 was obtained from feco in transmission. See table 2. For the 2-, 5-, and 20- 
mm spacers, the correct integers m/ were selected as described in section 5.1. 
Each set ol points, (m,-r-/,)X,-, lies along a curve parallel to the r 4 s curve. The 
constants subtracted from (77i/4-/y)\/ values for the 20-, 5-, andj2-mm spacers were 
401047-10 3 , 99552- 10 3 and 39592-10 3 A respectively. 
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Table 3 gives data for the calculation of wave- 
lengths 5570 and 4502 A, taking as the standard the 
line at 6056 A. Both method I and the "standard" 
method (eq 6) have been used. 

Table 3. Wavelength calculations from Fabry-Perot and feco 
data obtained with multilayer reflectors 

a. Calculation using eq (26), 20-mm spacer and feco. 



\s 


(OT.+/.) 


(m,-f-/«)X. 


T48(X S ) 


A 

6056.1620*... 


66276. 498 


A 

401,381.209 


.4 

333, 030 


M 


T48(X,) 


(«!«+/.) 


X, 


M 


A 

5570.32 
4502.38 ... 


A 

333, 650 
333, 070 


72057. 199 

89148.651 


.1 

*5570. 3224 
*4502.3816 


A 

5570. 3222 
4502. 3809 



b. Calculation 


using eq (6), 20- 


run and 5-mm spacers. 


M 


(«.'+//) 


(m x '+f x f ) 


X r 


A 

5570.32 

4502.38 


16493. 291 
16493.291 


17931.927 
22185. 150 


A 

*5570. 3215 
•4602. 3806 



*The wavelengths of the standard line at 6056 A and of the unknown lines 
are for the ambient conditions given in table 2. 
f Rough value. 
JFrom Bruce [31] (converted to ambient conditions). 



6.3. Precision of Measurements 

In the determination of wavelength of fringes the 
standard deviation for an individual reading was 
generally from 0.1 to 0.3 A. This resulted in a 
standard deviation for the values of r of from 5 to 
10 A, with a rare value as high as 30 A. With 
Ar=10 A as an example and A = 5000 A, the ratio 
2A</>/\ discussed earlier in connection with eq (9) 
becomes 2A<£/\=At/\= 0.002. This is to be com- 
pared with the usual error in exact fraction determi- 
nation, A/=0.01. 

The improvement in precision of wavelength 
measurement is not so dramatic, however, since this 
is limited by the errors in measuring the fractions 
f s and j x . For example, in the standard method, 
eq (6), we take data from the 20-mm and 5-mm 
spacers (table 2) on the known line and the unknown 
line, 5570 A. For A/ S =A/ X =0.01, eq (7) yields 
AX X =0.0022 A. 

To compute the error obtained in the method em- 
ploying feco, we take A/,= A/* = 0.01 and 2A<j>(\)/\ = 
Ar/X=0.002. Then eq (9) yields AX* = 0.0012 A. 

The precision of both methods can, of course, be 
improved by repeated readings. For example, 
Stanley and Meggers [34] report that in the measure- 
ment of iron lines with 20- or 25-mm spacers, the 
average of 5 or 6 readings had a probable error of 1 
part in 7 million. 

The question arises as to the best etalon spacing to 
be used for feco. Since the wavelength is multiplied 
by the order number to obtain the r curve, smaller 



spacings and the resulting low order numbers tend 
to give smaller errors. But smaller spacings give 
fringes widely spaced in the spectrum. If the phase 
shift has a nonlinear variation with wavelength, as 
with the broadband multilayer, it is desirable to have 
the experimental points fairly close together. In 
this case, the best spacing was found to be 10 to 20//, 
whereas for aluminum the smallest possible spacing 
was desirable. 

6.4. Differences Between Transmission and 
Reflection Data 

By placing a semireflecting mirror between the 
collimator and interferometer in figure 4 it is possible 
to project the light reflected from the interferometer 
onto the spectograph slit. With a small separation 
between the interferometer plates the resulting spec- 
trum consists of dark fringes crossing the bright 
continuum. See figure 5a. 

If the reflecting surfaces are nonabsorbing, the 
principle of conservation of energy requires that the 
reflected spectrum be the exact complement of the 
transmitted spectrum. Under these circumstances 
the wavelengths of the fringes should be the same in 
the two spectra. Then the r n curve plotted from 
the reflection data should also give the phase shift 
dispersion. 

If the reflecting films are partially absorbing, this 
conclusion does not necessarily hold. Holden [35] 
lias shown that there can be asymmetry in the in- 
tensity distribution of fringes obtained in reflection. 

Experiments were conducted to determine whether 
or not the r n curve obtained from reflection feco was 
the same as that obtained by transmission. Ex- 
posures were taken within a few minutes of each 
other, and the same areas of the reflecting films were 
used. As expected, the aluminum films showed a 
small but significant difference. As shown in figure 
8 the slopes are different, the maximum discrepancy 
between 4200 and 6500 A amounting to about 40 A. 
The reflection fringes showed a definite asymmetry, 
which made determination of their wavelengths 
more difficult. 

The two curves for the 15-laver multilayer showed 
a difference of about 150 A between 4800 and 6600 A. 
See figure 8. This result was unexpected because 
of the low absorption of dielectric multilayers. 
Transmittance and reflectance measurments on 
these multilayers indicated with the absorptance 
is no more than 2.5 percent at 4800 A and no more 
than 1 percent at 6000 A. However, no error is 
introduced if phase shift data obtained by trans- 
mission feco are used in connection with Fabry- 
Perot fringes obtained in transmission since both 
are governed by the interference eq (15). The 
reflection fringes are governed by a different 
interference condition [11, 35, 36]. 

6.5. Uniformity of Phase Shift Dispersion 

With multilayer films the reflectance and phase 
shift on reflection are generally more sensitive 
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Figure 8. Differences in data obtained by transmission, and 
reflection. 

The experimental points represent (n-\-i)\,\ where n Is the order number of the 
fringe at Xo, i is an Integer, and x, Is the wavelength of a fringe observed either by 

transmission or by reflect ion. In the ease of the l.vlayer films a relatively thick 
spacer (4^/x) was used and hinges were closely spaced. Only about one-third of 
the fringes were measured and plotted, in the ease of the aluminum films the 
vertical scale interval has been expanded to four times that of the mult ilayer in 
order to show more clearly the small but significant difference between the trans- 
mission and reflection data. Only the t ransmission /Vro should be used in connec- 
tion with Fabry-Perot fringes obtained by transmission. 

functions of the thickness of the individual layers 
than with metal reflectors. If the absolute value 
of the phase shift is not constant over the aperture, 
it has the same effect as nonflatness of the support 
plate. If the dispersion of phase shift is not constant , 
then the shape of the reflected wavefront will depend 
on the wavelength [37]. In particular, if two such 
mirrors are adjusted to he optically parallel when 
a green line is used, they may not be optically 
parallel in the red or blue. 

The phase shift dispersion was measured at six- 
points across the 15-layer films. Over a 40-mm 
distance there was a definite change in the shape 
of the curve, amounting to 300 A from the blue to 
the red end of tin* spectrum. For precision Fabry- 
Perot interferometry this variation is too large to 
be tolerated. For the experiments reported here 
a 10-mm-diam aperture was placed between the 
interferometer plates. The total variation in phase 
shift dispersion over this area was less than 100 A, 
which is the same order of magnitude as the departure 
from flatness of the support plates. 

Freshly deposited aluminum films showed very 
little variation in phase shift dispersion over a 
40-mm aperture, the maximum discrepancy being 
about 40 A over the region 4300 to 6600 A. 

These results indicate that multilayers for inter- 
ferometers in which several wavelengths are to be 
used should be checked for uniformity of phase 
shift dispersion. The use of feco provides a con- 
venient and precise method. 

7. Conclusions 

It has been found that measurements of reflection 
phase shift dispersion by means of fringes of equal 

543405—60 2 



chromatic order are useful at three stages in wave- 
length comparisons by Fabry-Perot interferometry. 
First, for multilayer reflectors it must be established 
that (he phase shift dispersion is sufficiently uni- 
form over the aperture. This can be done by obtain- 
ing the feco al several points across the aperture. 

Second, I he phase shift data can be combined 
with the method of exact fractions to yield the inte- 
gral order numbers. This is especially valuable in 
the case of multilayers where the dispersion of the 
phase 4 shift is high and where calculated phase shifts 
may he in error due to poorly known layer thicknesses. 

Third, the accuracy of measurements of unknown 
wavelengths can be improved somewhat. The basic 
reason is that for high reflection mirrors, wavelengths 
of feco can be measured to about ±0.002 order, 
whereas the error in determining fractions in Fabry- 
Perot patterns is about ±0.01 order. This differ- 
ence is due primarily to the fact that feco arc obtained 
from a very small area of the plate. Small depart- 
ures from flatness of the substrate 4 surface do not 
reduce the sharpness of the feco. On the other hand, 
the Fabry-Perot fringes result from interference over 
the entire 4 aperture, and variations in the surface 
reduce the fringe sharpness. In return for the gain 
in precision afforded by the feco the necessary pre- 
caution is that the variation in phase shift dispersion 
across the 4 aperture is sufficiently small. 
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9. Appendix 

9.1. Derivation of Equation (28) 

It is to be shown that for an arbitrary integer q 
and for integers, r h selected as described, the differ- 
ences (fj+jj) \j—r Q (A y ) are equal to 2(£— f). To 
both sides of eq (19) we add the quantity (q—n) X . 

2Xo=2*' + 2* (\,) + (g-ft) x . 
Adding (q—n) \ x to each side of eq (20) gives 
(2+1) X^^+20 (X0 + (2-^) Xi. 

Therefore the smooth curve connecting the points 
<ZV (2+1) Xi, etc., is the continuous function 

T g (\)=2f+2<l>(\) + (q-n)\. (30) 

Now from eq (22) we obtain the relation 

(r a +fa)K= (m a +fa)K+ {r a —m a )\ a 

=2*+20(X a ) + (r fl -m fl )X fl (31) 

The integers r a , r b , etc., are selected so as to 
satisfy eq (27). Substitution of eq (31) and the 
corresponding equation for X 6 into eq (27) yields 

(r a +f a )\ a -T q (\ a ) = 2t+24(\ a ) + (r a -7n a )\ a 

-[2f+20(X a ) + (2-n)X a ]=2f+20(X B ) + (r 6 -m»)X 6 

-[2f+2<KX & ) + (2-n)X & ] 

or (r a — m a — q+n)\ a = (r b —m b —q+n)\ b . (32) 



Since r, m, q, and n are integers, the only general i 
way in which eq (32) can be satisfied is r a —m a 
=r b —rn b = q—n, in which case 



(ri+/i)X i -T f (X i )=2(*-f) 



(33) 



where j stands for a, b, or c. Thus the difference 
between the etalon spacings can be determined, 
regardless of the phase shift dispersion and of the 
integer, q, which is arbitrarily selected. 



9.2. Required Accuracy of Theoretical Phase Shift 
Curve 

The phase shift curve must be known from theoret- 
ical calculations with a certain accuracy if the 
"true" order numbers are to be determined, as in 
section 5.1. We wish to show that the slope of the 
curve 20 versus X must be known to better than 
±0.5. 

With jeco the criterion for selection of the true 
order number, n, for the fringe at X is that the plot 
of nXo versus X , (w+l)Xi versus Xi, etc., gives a curve 
r n (X) which is parallel to 20(X), i.e., r n (X) = 2/' + 20 
(X). The curve which results for an incorrect 
integer, q, is given by eq (30). Now the smallest 
possible error is one unit, i.e., q—n-\-l. Then by 
eq (30), 

r n+1 (X) = 2/' + 20(X) + X. 

It is evident that at all values of X the slope'of the 
T n+ i(X) curve is greater by one than the slope of the 
r n (X) curve. Thus the possible r q curves are a 
family of curves which at a given value of X differ 
in slope by one. In order to make the correct 
selection of n, therefore, the slope of the 20 (X) curve 
should be known to better than ±0.5. 

Similarly for Fabry-Perot data, the criterion for 
selection of the order numbers m j (j=a,b 1 or c) is that 
when the values of (m ; +/ ; )X ; - are plotted against \ } 
the points fall along a curve, A, parallel to 20(X). 
If the integers, m h are each increased by one, the 
resulting points (m.,+ l+jQX ; - versus X ; fall along a 
curve having a slope greater by one than curve A. 
Therefore in order to make the correct selection of 
m h the slope of the 20 (X) curve should be known to 
better than ±0.5. 

This is not a severe restriction, amounting to a 
maximum permissible uncertainty in € = 27r0/X of 
about ±35° over the wavelength interval between 
5000 and 6000 A. Therefore in most cases of 
metallic reflectors and multilayers the theoretical 
phase shift will be known accurately enough to 
calculate the true order numbers. When it is not 
known with sufficient accuracy, method II can be 
used. 
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